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Phosphorus (P) adsorption and desorption characteristics of Min- 
nesota Basalt Lunar Simulant (MBLS) and Lunar Glass Simulant (LGS) 
were evaluated. Results of P interactions with lunar simulants indicated 
that mineral and glass components adsorbed between 50 and 70% of the 
applied P and that between 85 and 100% of the applied P was desorbed. 
The Extended Freundlich equation best described the adsorption data (r 2 
= 0.92), whereas the Raven/Hossner equation best described the desorp- 
tion data (r 2 = 0.97). Kinetic desorption results indicated that MBLS and 
LGS released most of their P within 15 h.The expanded Elovich equation 
fit the data best at shorter times while t/QoT equation had a better fit at 
longer times. These results indicate that P does not strongly adsorb to the 
two simulants and that any P that was adsorbed was readily desorbed in 
the presence of anion exchange resin. This work suggests that multiple 
small applications of P (10-20 mg P kg“ ! ) should be added to the simu- 
lants to ensure adequate solution P for plant uptake and efficient use of 
P fertilizer. 


C ROPS will be grown in lunar soil at an astro- 
naut tended-lunar base to produce food and 
oxygen, recycle carbon dioxide, and supply water 
(transpiration). Because plant essential nutrients 
will be a limiting component of a lunar Advanced 
Life Support System (ALSS), it is important to 
understand plant nutrient interactions with luTTar 
soil. High concentrations of calcium (1 13,640 mg 
kg' 1 ), magnesium (45,230 mg kg -1 ), and sulfur 
(300—1400 mg kg -1 ) are present in lunar soil 
(Hossner and Allen 1989). Phosphorus (P) has 
been reported to occur in lunar rocks and soils, 
with an average concentration of 800 mg P kg' 1 
(Hossner and Allen 1989). However, it is un- 
known how fertilizer P will interact with lunar 
soil when applied for plant growth. 

Secondary minerals such as phylosillicates, 
X-ray amorphous colloid material (e.g., allo- 
phane), iron and aluminum hydrous oxides, and 
calcium carbonate have been shown to be re- 
sponsible for P adsorption in terrestrial soils, 
which results in less P being available for plant 
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growth (Sample et al. 1980; Sanchez and Uehara 
1980). However, depending on the type of P 
bonding that occurs, P can be released (des- 
orbed) into solution in the presence of a sink 
(e.g., plant root). 

Quantity/ intensity parameters are utilized to 
describe P adsorprion/desorption relationships. 
Quantity (Q) refers to the amount of P adsorbed 
to the soil in the case of adsorption or the 
amount of P released from the soil in the case of 
desorption. Intensity (/) refen to the amount of 
P in solution. Phosphorus adsorption and des- 
orption plots are obtained by graphing O venus 
/. Adsorption or desorption buffering capacity 
can be obtained by plotting the derivative of the 
Q/7 plot venus /. Buffering capacity is defined as 
the amount of P that is added or removed that 
changes the solution concentration by one unit 
(Raven and Hossner 1994a). A more general de- 
finition of buffering capacity is the ability of a soil 
to maintain P solution concentration when P is 
added or removed from the soil system. 

Since lunar regolith is considered a national 
treasure and is only allocated in small quantities 
(mg), terrestrial analogs were utilized to under- 
stand how lunar soil interacts with plant nutri- 
ents. Research concerning plant interactions with 
lunar soil (Weete and Walkinshaw 1972; Walkin- 
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shaw etal. 1973; Walkinshaw and Johnson 1971; 
Baur et al 1974; Shay ct al. 1974) and chemical- 
mineralogical interactions in terrestrial analogs 
have been performed (Easterwood etal. 1991; 
Oglesby et al. 1993; Tan and Henninger 1993; 
Eick et al. 1996), but work concerning P interac- 
tions has not been done.The objective of this re- 
search was to determine the P adsorption, des- 
orption, and kinetic desorption characteristics of 
the Minnesota Basalt Lunar Simulant (MBLS) 
and the Lunar Glass Simulant (LGS). Quantity 
and intensity parameters were utilized to under- 
stand P interactions with lunar simulants. 

MATERIALS AND METHODS 
Lunar Simulants 

Minnesota Basalt Lunar Simulant (MBLS) 
was obtained from a 1 to 2-m-thick by 50-m- 
long sill found along a rock face at an abandoned 
quarry in Duluth, Minnesota (Weiblen and Gor- 
don 1988). This basalt is a suitable lunar soil sim- 
ulant because its bulk chemical and mineralogical 
composition is similar to Apollo 1 1 mare soil 
sample 10084 (Weiblen and Gordon 1988). 

Lunar glass simulant (LGS) was obtained from 
the Corning Glass Co. (Corning New York) and 
has a bulk chemical composition similar to the 
average of high titanium soil found at the Apollo 
1 1 landing site (C. Galindo, 1994, NASA, Johnson 
Space Center, personal communication). The 
Corning Glass Co. modified a procedure oudined 
by Minkin et al. (1976) to produce LGS. 

Both simulants were ground to have the same 
size distribution as Apollo sample 10084. How- 
ever, lunar simulants are not exact replicates of 
lunar soil. For example, they do not contain cer- 
tain components of lunar soil such as agglutinates. 
Agglutinates are particles made of lithic, mineral, 
and glass fragments melted together by heat gen- 
erated during meteorite impacts and are ex- 
tremely difficult and expensive to produce syn- 
thetically. Despite the lack of agglutinates, these 
simulants can be used to provide insight into P 
reactions in lunar soil. 

Characterization 

X-ray diffraction analysis of the simulants was 
performed using a Phillips diffractometer (Cu- 
Ka X-ray tube). Chemical analysis included pH 
(McLean 1982), electrical conductivity (Rhoades 
1982), exchangeable canons (Ca~ + , Mg" + , K + , and 
Na + ) (Thomas 1982). total P (Olsen and Som- 
mers 1982; lackson 1958), extractable P (1.4 M 
NH 4 OAc and 0.025 M EDTA at pH 4.2) (TAES 


1980), and oxalate (Sheldnck 1984) and sodium 
dithionite extractable iron (Mehra and Jackson 
1960). Microbiological analysis consisted of bac- 
terial (Zuberer 1994) and fungal (Parkinson 
1994) counts. Physical properties determined 
were size distribution (Jackson 1969;. wilting 
point, and field capacity (Cassel and Nielsen 
1986). 

Simulant Preparation 

Phosphorus was added as KH 2 PO solution 
to triplicate simulant samples m 5-mL 
minibeakers at rates of 0, 2.5. 5.0, 10, 1 5, 20, and 
25 mg P kg" 1 for the adsorption experiments. 
The 10 mg P kg -1 rate was utilized in the des- 
orption and kinetic desorption experiments. 
Anion exchange resin weights and d -sorption 
times were varied. Simulant weight for the MBLS 
adsorption experiments was 1 g. and 2 g of 
MBLS were used for the desorption and kinetic 
desorption experiments. The LGS experiments 
used 0.75 g of material because of a limited sup- 
ply. The P-treated samples were incubated in a 
glove box for 7 days for the adsorption and des- 
orption experiments, and the kinetic desorption 
experiments were incubated for 8 day ,. Samples 
were incubated at field capacity' (0.03 MPa mois- 
ture), 23 to 25°C, and 85 to 95% humidity. In an 
attempt to mimic the atmospheric conditions of 
the ALSS, the incubations were performed in 
C0 2 , 0 2 , and N 2 concentrations of 0.1, 21.75, 
and 78.15%, respectively (D.W. Ming, 1992, 
NASA, Johnson Space Center, persona! commu- 
nication). All P analyses used the stannous chlo- 
ride method (Jackson 1958). 

Phosphorus Adsorption 

Phosphorus adsorption experiments are 
modified slightly after methods described by 
Raven and Hossner (1994a). After incubation, 
the simulant was transferred into 150-mL 
beakers, with a 2:1 deionized water to simulant 
ratio. The samples were shaken for 1 h, and the 
solution was then filtered through a rmllipore fil- 
ter (0.45 pm) and subsequently analyzed for P. 
The amount of P adsorbed was determined by 
subtracting the amount of P in solution from the 
amount of native solution P plus added P. Native 
solution P was determined by analyzing the solu- 
tion with the procedure outlined above with the 
exception that no P was added to the simulant. 

The amount of P adsorbed (Q A ) was plotted 
versus solution P (I A ), and the Langmuir, Fre- 
undlich. Extended Freundlich, andTemkin equa- 
tions (Barrow 1978) were used to describe the 



Vol. 161 - No. 12 


Phosphorus Properties of Lunar Simulants 


875 


data. Data points associated with Q A = 0 were 
not used in the Langmuir and Freundlich equa- 
tion regression analysis because ln(0) fFreundlich) 
and I A /0 (Langmuir) are not real numbers. 

Phosphorus Desorption 

The Raven and Hossner (1993) method was 
utilized for the MBLS P desorption studies. 
Dowex 1-X8 saturated with HCO* 3 anion ex- 
change resin contained within nylon mesh (Nitex 
100% polyamide nylon fiber, mesh 3-210/36) 
was used in the desorption studies. The MBLS 
material was pre-equilibrated with 10 mg P kg 1 
simulant under the conditions mentioned. After 
equilibration, MBLS was transferred into 150-mL 
vials with a 2:1 water to simulant ratio and 
shaken on an orbital shaker for 1 h. After shak- 
ing, resin bags containing 0, 0.001, 0.0042, 
0.0075, 0.015, 0.025, 0.25, and 1 g of resin (dry 
weight) were added to triplicate simulant suspen- 
sions using 3 to 4 drops of chloroform to reduce 
microbiological activity. The suspensions were 
shaken for 72 h, and the resin bags were then re- 
moved and washed free of particles with deion- 
ized water. The resin bags were placed into 20 mL 
of 0.5 Af NaCl and shaken for 30 mm to desorb 
the P from the resin. The simulant suspensions 
were allowed to equilibrate for 1 h, centrifuged at 
3000 r.p.m. for 10 min, and filtered through a 
millipore filter (0.45 |lm). 

The amount of P desorbed from MBLS was 
determined from the resin extract, whereas the 
solution concentrations were determined as in 
the adsorption experiments. Native labile-? was 
determined by the procedure outlined above, but 
with zero P added and using 1 g of resin. The 
amount of P desorbed (Q D ), compared with so- 
lution P (1 D ), was fitted to the Raven and Hoss- 
ner (1993), Barrow (1979), and Brewster et al. 
(1975) equations (Table 2). When the Ip mea- 
surements were below the detection limit, 0.001 
mg P kg -1 was used in the regression analysis. 

Desorption Kinetics 

The same procedure used in the desorption 
experiments was used in the kinetic desorption 
experiments, except that 1-g resin weights were 
used and the time of desorption was varied from 
0.5 to 72 h. Triplicate samples for MBLS and du- 
plicates for LGS were used. 

Kinetic data were plotted as desorbed P 
(Q dt ) versus time (t).The first order, second 
order, two-constant (Elkhatib and Hern 1988), 
parabolic (Cooke 1966), and expanded Elovich 
(Raven and Hossner 1994b; Polyzopoulos et al. 


1986) equations were used to describe the kmedc 
data (Table 3). Also used was the t/Q L: r equation, 
which plots t/ Qdt versus t, yielding the linear 
form of t/QoT = d + &t, which converts to Qdt 
= t/(a + bt) (Table 3). Because t/Q D r was not a 
real number at Qdt = 0, a value ot 0.001 was 
used at 0 t and 0 mg P kg' 1 Q D t in the t/Q DT 
equation. The parameters (a,b) for the t/Q D T 
equation came from the linear form Once the 
parameters were obtained, they were used in the 
Qdt = t/(a +bt) equation, which was plotted and 
compared with the original data. 

RESULTS 

Characterization 

Mineralogical analysis (XRD) indicated that 
MBLS contained anorthite, augite, and llmemte. 
This is in agreement with other mineralogical 
analyses of MBLS (Tan and Hennmger 1993; 


TABLE 1 


Characterization data for Minnesota basdt lunar 
simulant (MBLS) and lunar glass simulant (LGS) 



MBLS 

LGS 

Size Fractions (mm) 

% WT. 

Sand>50 

54.3 

54.1 

Silt 2-50 

44.2 

44.2 

Clay<2 

1.5 

1.7 

Moisture (MPa) 

% wt. 

0.03 

8.5 

9.4 

1.5 

0.9 

1.5 


Exchangeable Cation cmol kg ! 


Ca 

0.53 

0.47 

Mg 

0.21 

0.17 

K 

0.06 

0.04 

Na 

0.14 

0.02 

pH 

9.30 

9.75 



S 1 

E.C. 

0.007 

0.008 



mg kg 

Oxalate Fe 

12115 

1684 

DCB Fe 

3715 

359 

Total P 

140 

417 

Extractable P 

28 

t 

Resin Desorable P 

1.18 

0.78 

Microbes 


counts ki ”* 

Bacteria 

1X10 5 

1X10 5 

Fungi 

1X10 5 

1X10 5 


fNot determined. 
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Gordon 1991). Because of its glassy nature, the 
LGS material had no prominent X-ray peaks. 

The simulants were ground to the same size 
distribution as Apollo 11 sample 10084. The size 
fraction data (Table 1) indicate that the simulants 
have a sandy loam texture according to the USDA 
textural triangle (Singer and Munns 1987). Be- 
cause of the sandy nature of these simulants, field 
capacity (0.03 MPa) values were found to be low 
for both simulants (Table 1). 

The elevated pH values found in MBLS and 
LGS (Table 1) can be explained by hydrolysis re- 
actions of the primary minerals in MBLS (Stevens 
and Carron 1948; Stumm and Morgan 1981; 
Bohn et al. 1985) and alkali and alkaline cations in 
LGS. The hydrolysis of the MBLS mineral anor- 
thite can be illustrated as follows (Stumm and 
Morgan 1981): 

CaAl 2 Si 2 0 8(l) + 3H 2 0 — > 

Ca 2+ + 20H“ + Al 2 Si 2 0 5 (0H) 4l 

It should be noted that this model is being used 
solely to depict the increase in pH associated with 
the weathering of anorthite. 

The term abrasion pH is defined as the pH ot 
a mineral solution after that mineral has been 
ground in water (Stevens and Carron 1948). As 
the primary minerals are ground in water, hy- 
drolysis reactions occur resulting in elevated pH 
values. Because the preparation of MBLS re- 
quired grinding, similar effects associated with 
abrasion pH determination appear to have caused 
high pH values in MBLS. 

These simulants were unweathered scrthey 
were not highly charged and, therefore, possessed 
low amounts of exchangeable cations (Table 1). 
Broken edges may have provided enough edge 
charge for limited cation exchange.The electrical 
conductivity (EC) of both simulants was low 
(Table 1), indicating that these simulants would 
not pose a salinity problem. 

Dithionite-citrate-bicarbonate (DCB) ex- 
tractable Fe values were lower than the oxalate 
extractable Fe (Table l).This is contrary to ex- 
tractable Fe values for most terrestrial soils. The 
DCB method is more suitable for removing oxi- 
dized Fe, and the oxalate method is more suitable 
for removing reduced and poorly crystalline Fe 
(Wang et al. 1991). Because the simulants were 
composed mostly of reduced Fe (MBLS) and 
poorlv crystalline Fe (LGS) (Galindo et al. 1989), 
the DCB method was not expected to extract 
more Fe than the oxalate method. 

Total P concentrations in the surface 8 inches 
of terrestrial soils have been found to range from 


38 to 956 mg P/kg, with an average of 382 mg 
P/kg (Tisdale et al. 1985). The MBLS material 
fell below the average, whereas the LGS material 
exceeded the average (Table 1). Extractable P 
(Table 1) for MBLS was in the medium range. 
Resin-extractable P (Table 1) was higher in 
MBLS than in LGS. Because total P ww higher in 
LGS than in MBLS, it appeared that LGS retained 
P more strongly than did MBLS. 

Because the simulants lacked orgamc matter, 
microbial counts were low (Table 1) when com- 
pared with most soils (10 8 - 10 9 bactena kg’ 1 , 10 8 
- 10 n fungi kg’ 1 ) (Alexander 1977). 

PAdsoqjtion 

The data indicate that added P was not ad- 
sorbed strongly by the simulants (Fig. I ). From 80 
to 100% of the added P was adsorbed at applica- 
tion rates below 5 mg P kg' 1 . Only 50 to 70% of 
the applied P of the applied P was adsorbed at 
the 25 mg kg* 1 treatment level. Most terrestrial 
soils would have continued to adsorb higher per- 
centages of P beyond the 25 mg kg treatment 
level. Texas soils, ranging from acid to slightly 
basic (Raven 1992), adsorbed 78 to 95% of P 
added at a rate of 25 mg P kg" 1 .Volcanic ash soils 
from Guatemala retained 92 to 97% of the 929 
mg P kg -1 applied P (Prasad and Motto 1982). 
Subtropical acid soils from Australia adsorbed 
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Fig. 2. Extended Freundlich plots of Minnesota basalt 
lunar simulant (MBLS) and lunar glass simulant (LGS). 
Error bars represent standard deviation. Error bars on 
some points are to small to be observed 


Fig. 3. A two-region Langmuir plot of U/Q A versus l A for 
Minnesota basalt lunar simulant (MBLS) and lunar glass 
simulant (LGS). Error bars represent standard devia* 
tion. 


100% of added P (25 mg P kg~ l) . When com- 
pared with other soils, it is evident that P was not 
strongly adsorbed to the simulants. However, 
sandy soils do not adsorb high amounts of P be- 
cause of the lack of P-adsorbing minerals (sec- 
ondary’ minerals). Secondary minerals (e^., Fe 
and Al-hydrous oxides), which are largely respon- 
sible for P adsorption, were not detected in the 
simulants; therefore, P adsorption was not exten- 
sive. 

All adsorption equations described P adsorp- 
tion to the simulants (Table 2). However, the Ex- 
tended Freundlich best described the adsorption 
data for both simulants (Fig. 2) because of the 
hieher r value. This indicates that the energy of 
adsorption decreased exponentially with increas- 
ing P adsorption (Bache and Williams 1971). The 
Extended Freundlich and the Freundlich equa- 
tions both described the data equally well for 
MBLS, but the Extended Freundlich equation is 
favored because it accounts for native adsorbed P 
in the c term. The maximum adsorption capaci- 
ties for MBLS and LGS as determined by the 
Langmuir equation (a term; Table 2) (Barrow 
1978) were 17.2 mg P kg~ 1 and 19.2 mg P kg \ 
respectively. When the P adsorption data were 
plotted in the Langmuir form (1 A / Qa vs Ia), two- 


region Langmuir plots were obtained (Fig. 3). It 
has been proposed that the two-region Langmuir 
plot indicates adsorption (region 1) and then pre- 
cipitation (region 2) (Griffin and Junnak 1978). 
Therefore, it is possible that P adsorption oc- 
curred up to about 0.40 mg P/L in MBLS and 
0.3 mg P/L in LGS. Beyond these L values, the 
concentration of P in solution was high enough 
to cause precipitation of calcium phosphate. 

The adsorption buffering capacity of both 
simulants was considered low when compared 
with selected Texas soils. Buffering capacities ot 
Texas soils ranged from 37 to 252 L kg 1 at 0.1 
mg P/L (Raven 1992), whereas MBLS and LGS 
had values of 20 L kg 1 and 24 L kg \ respec- 
tively (Fig. 4). A low adsorption buffering capac- 
ity indicates that these simulants have a limited 
ability to adsorb P and that adding P will result in 
elevated solution P concentrations. 

Desorption Study 

Desorption data for MBLS shoved that 88% 
of total P (total P = native soluble P + applied P) 
was recovered with 1 g resin (Fig. 5) At the low- 
est resin treatment (0.025 g), 67% of the P was re- 
covered from the MBLS material. Selected sandy 
soils from Niger (Wendt et al. 1993) released 28 
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TABLE 2 


Equation parameters and statistical data describing P adsorption onto 
Minnesota basalt humor simulant (MBLS) and lunar glass simulant (LGS) 



MBLS 

LGS 


Parameters 

Statistics 

Parameters 

Statistics 

Extended 

a = 13.0 

r 2 = 0.97 

<1 = 15.2 

= 0.92 

Freundlich* 

6 = 0.5 


b = 0.4 


In(Q A +c)=liw+Mn (I/O 

C = 1.2 


c = 0.8 


Freundlich 

a = 11.5 


a = 13.9 


ln(Q A )-lna+Mn(lA) 

6 = 0.6 

r 2 = 0.97 

fe = 0.4 

r = 0.90 

Langmuir* 

a= 17.2 


a = 19.2 


Ia/Qa=1/<* b+bl/ab 

6 = 2.0 

r 2 = 0.89 

6 = 3.3 

r~ - 0.78 

Temkin 

a = 2.9 


a = 2.3 


Qa=<> In (b 1/J 

b = 46.0 

r 2 = 0.90 

b = 264.0 

r 2 = 0.74 


Prob>F=0.0001 
t c= native adsorbed P 

5 <i=max adsorbed P, ^binding energy constant 


to 81% of applied P to resin from 10 mg P kg" 1 
treatments. It should be noted that P recovery by 
Wendt et al. (1993) consider only the percent of 
applied P desorbed. Raven and Hossner (1993) 
obtained 67 to 79% desorption of applied and na- 
tive desorbed P in the range of 0 to 100 mg P 
kg' 1 . Calculations in this paper consider native 
desorbed P as well as added P when determining 
the percent P desorbed. Desorption data sug- 
gested that MBLS tends to release a higher per- 



Fig. 4. Adsorption buffering capacity (6QJ6\ A ) plots for 
Minnesota basalt lunar simulant (MBLS) and lunar glass 
simulant (LGS). 


centage of applied P than do terrestrial soils. 

Native P desorbed was determined to be 1.18 
mg P kg _l for MBLS. Because MBLS did not 
desorb more than 11 mg P kg' 1 at the 10 mg P 
kg” 1 treatment level, it was evident chat not all 
adsorbed P was recovered (Fig. 5). Barrow and 
Shaw (1977) reported that anion-exchange resin 
cannot extract all desorbable P from soil. 

The Brewster et al. (1975), Barrow (1978), 
and Raven and Hossner (1993) equations were 



Fig. 5. Ravcn/Hossner and Barrow equations fitted to 
Minnesota basalt lunar simulant (MBLS) P desorption 
data. Error bars represent standard deviation. 
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TABLE 3 


Equation parameters and statistical data describing P desorption from Minnesota basak 
lunar simulant (MBLS) for the Raven/Hossner, Barrow, and Brewster equations 


Equation 

Parameters 

Statistics 

Raven/Hossner 

a = 84.3 


Qo = dP""’ + Wn(In + 1)+r 

b = -44.7 

R ; = 0.97 


r = -80.7 


Barrow 

*1 = to.o 


Qd = ^ + M d c 

b = 31.1 

ii 

p 

‘o 


f = 0.5 


Brewster 

d= -1.5 


Q D = a + Wn(I D ) 

T 

II 

-45 

r - j 

QO 

d 

il 


Prob > F- 0.0001 


utilized to describe the MBLS desorption data. 
The Barrow and the Raven and Hossner equa- 
tions had the higher r 2 values (Fig. 5; Table 3). 
The Raven and Hossner equation would better 
describe desorption of P because it does not 
allow solution P to be reduced to 0 mg P/L. 

The Raven and Hossner equation was used 
to determine the MBLS P desorption buffering 
capacity. By taking the derivative of the Raven 
and Hossner equation (dQ D /dI D = -0.01 alo'* 1 01 
+ b/(I D +1)) and plotting it versus I D , a desorp- 
tion buffering curve was developed (Fig. 6). 

Buffering capacities of 480 L kg -1 at 0.01 mg 
P L _1 and 234 L kg -1 at 0.10 mg P L“ l were re- 
ported for selected Texas soils (Raven 1992). 
These buffering capacity values are 348 (0.01 mg 
P L~ ! ) and 185 L kg -1 (0.10 mg P L” 1 ) higher 
than buffering capacity values of MBLS at those 
corresponding solution P values. The lowffesorp- 
tion buffering capacity indicates that MBLS has a 
much lower ability to maintain solution P con- 
centration as P was removed from solution com- 
pared with terrestrial soils.The MBLS desorption 
buffering curve indicates that MBLS easily re- 
leased adsorbed P to the resin. 

Kinetic Desorption Study 

Most of the P was desorbed from both simu- 
lants within 15 h, similar to most terrestrial soils 
(Figs. 7 and 8). The MBLS material desorbed all 
of its P within 15 h, whereas LGS continued to 
desorb slightly more P with time. Most soils 
rapidly release P within 20 h and then continue 
to release P slowly over time (e.g., 90 h) (Raven 
and Hossner 1994b; Pavlatou and Polyzopoulos 
1988; Evans and Jurinak 1976). The high total P 
(417 mg P kg' 1 ) in LGS apparendv allowed con- 
tinual desorption of P beyond 15 h. 

The solution P values for both simulants were 
reduced to undetectable levels at 5 and 24 h for 


MBLS and LGS, respectively. A slightly higher 
adsorption capacity and higher total P content al- 
lowed LGS to maintain solution P tor a longer 
time than did MBLS. 

The expanded Elovich and t/Q^j equation, 
first and second order, parabolic, and two con- 
stant equations were fitted to the MBLS and LGS 
kinetic data. The expanded Elovich and the 
*/Qdt equations were found to have the higher r 
values and, therefore, best describe the kinetic 
data for both simulants (Table 4 and 5). The 
t/ Qdt equation possessed higher r J values than 
the expanded Elovich for both simulants. The ex- 
panded Elovich equation described the data best 



Fig. 6, Phosphorus desorption buffering capacity 
(dCfo/dlo) calculated from the RavervHossner desorp- 
tion equation versus P remaining in soiubon (l D ) for Min- 
nesota basalt lunar simulant. 
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Fig. 7. Kinetic P desorption data showing the Expanded 
Elovich and t/Qor plots for the lunar glass simulant (LGS) 
with amount of P desort>ed (Qot) as a function of time 
(t). Error bars represent standard deviation. Error bars 
on some points are too small to be observed. 


Fig. 8. Kinetic P desorption data showing tne Expanded 
Elovich and t/Oo T plots for the Minnesotd basalt lunar 
simulant (MBLS) with amount of P desorbed (Qcr) as a 
function of time (t). Error bars represent standard de- 
viation. 


at shorter times (<24 h LGS; <30 h MBLS), 
whereas t/Q DT was better at longer times (>24 h 
LGS; >30 h MBLS) (Figs. 7 and 8). 

DISCUSSION AND SUMMARY 

The adsorption data indicated that the simu- 
lants did not adsorb applied P readily. The Ex- 
tended Freundlich equation best described the P 
adsorption data, indicating that the energy of ad- 
sorption decreased exponentially with concentra- 
non. Interpretation of the Langmuir plots sug- 
gests that P was adsorbed onto the simulants up 
to about 0.40 and 0.30 mg P L 1 in MBLS and 


LGS, respectively, and that P began precipitating 
beyond these solution P concentrations. The ad- 
sorption buffering capacity data show that both 
simulants possess low adsorption buffering capac- 
ities, indicating that added P will not be adsorbed 
strongly and will result in elevated solution P 
concentrations. These data agrees well with P ad- 
sorption data of sandy soils, where low additions 
of P are required to obtain optimal solunon P 
concentrations (Olsen andWatanabe 1970: Wendt 
et al. 1993) 

The MBLS desorption data were best de- 
scribed by the Raven and Hossner equation.The 


TABLE 4 

Parameters and statistical data for Minnesota basalt lunar simulant (MBLS) and lunar glass simulant (LGS) kinenc 
desorption equations. Coefficient of determination (r 2 ) for the t/Q D T equation is from its linear f ?rm 


Equation 

MBLS 

LGS 

Parameters 

Sumacs 

Parameters 

Suosocs 

Expanded 

a = 0.7383 


a = 0.8956 


Elovich 

b = 0.0001 

r 2 = 0.94 

b = 0.0001 

r 2 = 0.96 

Q nx = a In(t<-fc)+c 

f = 6.8964 


c = 7.7796 



d = 1.0000 


d = 1.0000 


t/ Q dt 

a = 0.1716 


a = 0.3185 


*/Qdt = ( d + ^t) 

b = 0.1022 

r = 0.98 

b = 0.0778 

r 2 = 0.97 


Prob > F = 0.0001 
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TABLE 5 


Statistical data for Minnesota basalt lunar simulant (MBLS) and lunar glass simulant (LGS) kinetic desorption equations 


Equation 


MBLS 



LGS 


r 2 


Prob > F 

r z 


Prob > F 

First order 

0.31 


0.0018 

0.63 


0.0019 

Second order 

0.00 


0.7965 

0.01 


0.7115 

Parabolic 

0.48 


0.0001 

0.63 


0.0007 

Two constant 

0.72 


0.0002 

0.83 


0.0001 


anion exchange resin desorbed a high percentage 
of P from the MBLS, indicating that P adsorption 
to the simulants is mosdy reversible. The desorp- 
tion buffering capacity indicated that MBLS des- 
orbed most of its adsorbed P to the solution for 
uptake by the resin. 

Kinetic data show that MBLS and LGS re- 
leased most of the adsorbed P within 15 h. The 
LGS material desorbed more P than was added, 
indicating that some soluble P was being supplied 
by the crystal phase. 

X-ray photoelectron studies have shown that 
iron (Fe) and aluminum (Al) hydroxides have 
constituted concentrations up to 4.2% and 12%, 
respectively, in a wide range of basalts (White and 
Hochella 1992). Because surface adsorption of P 
onto mineral surfaces is apparendy attributable to 
iron and aluminum hydroxides, the low concen- 
trations of surface Fe and Al hydroxides would 
not allow large scale P adsorption. 

Unweathered volcanic ash soils (noncrys- 
talline colloids) with high Si solubility have been 
shown to have low P adsorption capacity' (Fox 
1974). As the volcanic soils weathered and the Si 
solubility decreased, P adsorption increased in 
these volcanic soils (Fox 1974). Similar to un- 
weathered volcanic soils, the low P adsorption ca- 
pacity of LGS may possibly be a consequence of 
soluble Si, which forms silicate species that block 
P adsorption sites (Clark and McBride 1984;Eick 
et al. 1996). A low concentration of Fe and Al at 
the surface may also contribute to the limited P 
adsorption capacity of LGS. 

The data indicate that P will interact with the 
mineral and glass fraction of lunar soil in a similar 
manner if P is added to lunar soil in an ALSS en- 
vironment. There are other components of lunar 
soil (e.g., agglutinates) that are not represented in 
the simulants, and their reactions with P are un- 
known. However, by attempting to simulate a 
component of lunar soil, we can begin to under- 
stand how P interactions with lunar soil proceed. 
The data suggest that obtaining the most efficient 
P fertilization of lunar soil would require a small P 


application (e.g., 10-20 mg P kg -1 ) before plant- 
ing, followed by similar treatments applied during 
crop growth. This fertilization scheme is based on 
wheat growth experiments (data not shown) and 
will vary according to the crop to be grown. This 
recommendation agrees well with the optimal P 
application of 10.2 mg P kg' 1 required for millet 
growth in sandy soils of Niger (Wendt et al. 
1993). Because the simulants had a low P adsorp- 
tion capacity, one large P treatment (e.g., 100 mg 
P kg -1 ) to lunar soil may put more P in solution 
than is required for plant growth, thereby pro- 
moting inefficient use of P fertilizer and nutrient 
imbalances (e.g., Fe, Mn, Zn) in the plant. As 
crops are grown in the lunar soil, the formation of 
secondary minerals (e.g., Fe-hydroxides, clays) 
will result (Easterwood et al. 1991), and the P ad- 
sorption and desorption characteristics of lunar 
soil will change over time. Therefore, understand- 
ing P interactions in lunar soil will be an evolving 
challenge for agronomists. 
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